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Abstract 

We address the problem that the cross section for the collisions of un- 
stable particles diverges, if calculated by standard methods. This problem 
is considered for beams much smaller than the decay length of the unstable 
particle, much larger than the decay length and finally also for pancake- 
shaped beams. We find that in all cases this problem can be solved by 
taking into account the production/propagation of the unstable particle 
and/or the width of the incoming wave packets in momentum space. 



1 Introduction 



When one applies the Fcynman Rules and the Golden Rule to a collision of 
unstable particles, the cross section turns out to diverge. We summarize how 
this happens. The divergence occurs for instance in the graph 




Mb: 



(1) 




The lower half of this diagram looks like the decay of the muon. The consequence 
is that the momentum k may be on its mass shell. After all, that is what one gets 
from the decay of a muon: a muon neutrino on its mass shell. The propagator 
of this muon neutrino contains a factor 

that gets squared and integrated over as the Golden Rule tells us to do. Because 
the neutrino that contributes the just given factor can be on its mass shell, we 
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find that in lowest order in e the result will go as e~^. This is a divergent 
quantity since e is taken to be infinitesimal. This problem is of relevance for 
muon colliders as was already noted in [1]. 

The question that naturally poses itself is if the cross section really diverges, 
and if so, how this can be regularized, and if not, how we can compute it. In 
general we can say that there are two possibilities to solve the question of the 
divergent cross section. The first is that one takes into account that unstable 
particles cannot really be in or out states. This can be done by considering 
Feynman graphs that include the production process of the unstable particle. 
The second possibility is that one takes into account that the incoming wave 
packets are not really sharp momentum states but that we always have inter- 
ference between states with the same total momentum but with some of the 
momentum moved from one incoming particle to the other. This is of impor- 
tance if the peak structure of the matrix element is sharper than the size of the 
incoming wave packet in momentum space. We indeed have this in the case of 
the unstable particle, since the propagator that causes the divergence has no 
particular size attribute, except for the e, which is taken to be infinitesimal. 

This problem has mainly been discussed for the case that would be applicable 
to, for instance, muon colliders. Cf., [3], [6], [2]. We will mostly address a case of 
academic interest that has been considered alongside the realistic case, namely 
that of very wide beams (i.e., much larger than the decay length of the muon). 
We will try to apply the methods applicable to the realistic case also in the case 
of infinitely wide beams. For this we will need some modifications to the realistic 
case. This case is considered in the next section. We reconsider it because we 
use different, more covariant notations, best introduced in a somewhat more 
familiar setting, and also because we have some points to add. 

Finally we consider the case of pancake-shaped beams (i.e., much larger than 
the decay length in the transversal directions and much smaller than that in the 
longitudinal direction). 

2 The Realistic Case 

Here we discuss the solution as would be applicable in muon colliders, that 
is, for beams of which the size is much smaller than the decay length of the 
unstable particle. In this case the propagator of the unstable particle confines 
its momentum to the mass-shell and because production and collision can be 
chosen to be well-separated in space, we need not worry about the production 
process. There is, however, still the point that the matrix element is peaked 
sharply enough to notice interference between states where some momentum 
is shifted from one incoming particle to the other. This case has been solved 
in [2], however we again present this here using notations that are manifestly 
covariant. 

The quantum distribution function n{p,r), introduced in [3], is used to de- 
scribe the particles. It is defined to be given by 

n{p, r) = ^ 1 d'A 6{p ■ A) <p{p + i A) <p{p - i A)* e"'^-. (3) 

The S{p- A) is used to confine the particle to its mass shell. This implies that we 
use the approximation that the components of A are much smaller than the ones 
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of p. The quantum distribution function contains as much information about 
the state of a particle as a density matrix. From it the probabihty densities in 
momentum and position space can be found. They are given by 

P{p) = — d^rn{p,r); p{r)= dp—n{p,r). (4) 
m J J m 

Both are (in the approximation that they are sharply peaked in momentum 

space) zcroth components of four-vectors. This should, of course, be the case 
with densities. The probability measure that belongs to them is respectively dp 
(which is by definition equal to (2^^^^) and d^r. This looks pretty much as if 
n{p, r) were a joint probability for r and p, but of course such a thing cannot 
really exist in quantum mechanics and actually n{p,r) does not need to be a 
real function, which discards it as a probability density. 
We extend the definition of the luminosity to 

dL{pi,p2,p) = dpidp2 / d'^rni{pi,r)n2{p2,r + p). (5) 

mim2 J 

This reduces to the normal definition for dL by setting p = 0. The GD that 
occurs here stands for "Gramm determinant" and is equal to plp2 — (pi ■ P2)^- 
Also the definition of the cross section is extended, namely to 

d(T{pi,P2, A) = --yl==(27r)^5(pi +P2-q)dqi--- dqn 

4V — tjrU (6) 

M{pi + iA,p2 - ^A,q)M{pi - iA,p2 + |A,g)*. 

This gives the well-known definition for the cross section if we take A = 0. The 
exact expression for the number of events W (see for instance [5]) can, using 
these definitions, be written as 

W = JdL{puP2,p)-^d^pe-'^-Pda{pi,p2,A). (7) 

Normally one assumes that da does not vary much with A and for that reason 
it is safe to put A = 0. Then the integrals over A and p become trivial and 
the familiar result that "number of events is cross section times himinosity" is 
obtained. In our case, however, a has a pole at A = so this approximation 
cannot be made. 

The reasoning about how to get a finite cross section is essentially identical 
to what is given in [2], so we will not present it here. In the end the result is 

a = aTv j daredj^^S{k'^ -m^), (8) 

where the "red" in ddred stands for "reduced" meaning that the factors l/{k'^± 
ie) that cause the divergence have to be left out. a is the transverse size of the 
beam. We take this to be given by a = ^/tt a = ^/tt ■ 10 /im, where (for roimd 
beams) a is the standard deviation in the position of the particles in the beam in 
either direction perpendicular to the beam. This was also done in [2]. Because 
this part of the cross section is proportional to this size, the effect is called 
"linear beam size effect". If there are contributions to the cross section from 
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Figure 1: da/dt (fb/GeV^) vs. t (GeV^) for ^ = 100 GeV 



parts of phase space away from the singularity they should be added separately 
to the cross sections. These parts do not depend on the beam size. Below we 
will see that such a separation arises rather naturally from the shape of the 
graph of da/dk^. The part of the cross section that has to do with regions of 
phase space away from the singularity will be called the "regular cross section" , 
while the part that comes from regions near the singularity (or singularities) 
will be called the "semi-singular cross section" . 

As in [2], we consider the process /z~ + /z+ — > + e~ + Pg- For da/dt at 
100 GcV we find the graph in figure 1. The infinite spike at i = is caused by 
the instability. Figure 2 shows a detail of the same graph. Now the singularity 
is prominently present. The reader can compare the numbers on the axes of 
both graphs to get an idea how narrow the singularity actually is. To be able 
to calculate the regular part of the cross section without having to worry about 
the singularity, the cut t < —m^ is used (as was done in [2]). The two graphs 
just shown, justify this cut (i.e., ~ 0.01 GeV^ which is in the neighbourhood 
of the minimum of da/dt). 

The total cross section for pseudo-singular and regular cross sections are 
plotted in figure 3. For the regular cross section we find the same graph as in [2], 
but for the semi-singular cross section our graph is about a factor 1.7 higher. 
The reason for this difference appears to be twofold. In the first place, we did 
our calculations from standard model coupling constants, while [2] expresses 
the cross section in other decay widths and cross sections. If we take this into 
account, the factor 1.7 becomes a factor 2. This factor 2 is then due to an error 
in equation 46 in [2]. This equation should have an extra factor 2 on the right 
hand side. 

Consequently we find that the semi-singular cross section dominates up to 
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Figure 2: da/dt (fb/GeV^) vs. t (GeV^) for = 100 GeV around t = 0. This 
is a detail of the previous plot. 
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Figure 3: a (fb) vs. ^/s (GeV). The solid line is the regular cross section and 
the dashed one is the pseudo-singular cross section, a = ^ ■ 10 jim. 
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about 113 GeV. [2] has 105 GcV. Furthermore it should be noted that at about 
90 GeV the semi-singular and regular cross sections are about the same order of 
magnitude because the regular one has a peak there, caused by the Z particle. 
For a bit above threshold the semi-singular cross section dominates strongly. 
At about, say, ~ 85 GeV one can safely forget about the regular part. Above 
\fs ~ 150 GeV the pseudo-singular cross section does not play any role anymore. 
This is the case up to arbitrary high energies because asymptotically the semi- 
singular cross section goes down as XjiySs/s), while the regular one goes down 
as 1/s. 

To calculate these cross sections, six diagrams involving 7, and 
as fundamental bosons were added. The algebra necessary was done by the 
C-I-+ computer algebra library GiNaC which is described in [9]. After that 
the integrations were carried out by adaptive Simpson integration. Unstable 
intermediate particles were given propagators using the iMV prescription. This, 
of course, raises the issue of gauge invariance. It was checked that for high 
energy the cross section goes down as 1/s, but there may still be imnoticed 
terms suppressed by a factor F^/M^ that did not show up in these calculations. 
The authors hope to address these issues in a subsequent paper. 

In practice one docs not need the cut-off t < —mf^. This is because this 
cut-off is already implied by the cut-offs imposed by measurability. If we take 
\/i = 100 GeV and demand that the energy of the outgoing electron of the 
pseudo-singular process is at least 1 GeV and that the angle under which this 
electron appears is at least 2° we have made it impossible for the muon neutrino 
to be on its mass shell (actually we then find that t indeed is always negative 
and has an upper limit of about — 5to^). 

However, if one wants to include the linear beam size effect in the normal 
Monte Carlo integration procedure, one can do this by doing the replacement 



fc2+ie k^ + i\kjL\/a' ^ ' 

This, in spite of its ad-hoc appearance, gives precisely the correct answer. It 

effectively introduces a decay to the muon neutrino to account for the fact that 
it can no longer collide after it has left the beam. This way of handling the 
divergence supposes that number of particles in the beam goes to zero expo- 
nentially as wc leave the beam (i.e.. the iMT prescription is about exponential 
decay), which is not terribly realistic, but in most cases this will not matter. 



3 Infinitely Wide Beams 

3.1 The Wrong Way 

The method of calculation introduced in [4] (to be called Ginzburgs method in 
the rest of this section) starts by observing that the propagator of an unstable 
particle is given by {p^ — + imr)~^ . We then "conclude" that the mass of the 
unstable particle has acquired an imaginary part. The four momentum-squared 
of this particle should be complex too. In its rest frame it is given by 

pi = (TO-iF/2,0). (10) 
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Wc now write for the "new" value of k"^ , that is the value that takes the complex 
momentum components into account, 

^new = (Pi - 1? =rr?- + g2 _ 2(m - ir)g«, (11) 

where q = qi + q2 a-s drawn in equation 1. For some reason we take the same 
values for the components of q but only change the components of pi and k. 
Normally (i.e., without complex momentum components), the value of g° is 
given by 

After substituting this, and neglecting the small quantity proportional to fc^F, 
we obtain 

fcncw = fcold-«^("^^-9^)- (13) 

This complex value then replaces the one given in equation 2, and a finite result 
is obtained. 

The problem with all this is, of course, that it is not terribly difficult to 
think of a process where one of the outgoing/incoming particles involved is un- 
stable and then an incoming complex momentum component has, by momen- 
tum conservation, no place to go. Significantly, we never hear about momentum 
conservation for the other vertex in the diagram. 

In [6] it is shown that the result obtained by using the here described method 
is exactly what one would expect for a muon that decays in a medium of anti- 
muons. This result is 

sine d9d(j) , , ^, 

cr = {1 - cos6)w[lo)(Lj — (7^i^^w{s^f^). (14) 

The factor w{u)) du) sva.6 dO dcjj/A-K is the probability measure of finding a muon 
neutrino with a certain momentum. This calculation is done in the rest frame 
of the decaying muon. However, three remarks are in order here^ 

1. To obtain this result, the definition of the quantity cross section has to be 
modified; 

2. The same modification of definition can be aplied to our result (i.e., equa- 
tion 26, to be derived in the next section) and the result will be the same; 

3. It is a bit of a coincidence that the modified cross section of [6] turns out to 
have the same value as Ginzburgs method 

Let us disciiss these points in this order. 

Firstly (1), the number of events W is related to the cross section via 



W 



where is the four- volume in which the beams overlap and Ji_2 denotes the 
four-flux of the two beams. The space integral over Ji is the number-of-particles 
four- vector A''^ defined by 

TV" = TVu", (16) 
with N the number of particles. We thus have 



dW 

~dr 



^{N.-J^r^N^J^a^,. (17) 



^An email discussion with V.G. Serbo was helpful! to get these points completely clear. 
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On the other hand, we expect to be able to calculate the number of events from 
considering collisions between muon neutrino's and muons, taking the momen- 
tum distribution of the muon neutrino's into account. Doing so we find 

^ = y ^oHda; ^^"^^^^^^ ^(jV.(a;) • J^f - iV.MVfa.^sH) (18) 

After specializing to the rest frame of the /i~ we see that the division of the two 
flux factors gives the 1 — cos 61 (For this one has to assume that k"^ = ~ 0.) 
so that we indeed find equation 14 back. The modification to the definition 
of the cross section is that one equates these two dW/dfs. The first dt refers 
to the time the ^~ track spends in the /i+ cloud while the second dt refers to 
the time the decay product spends in this cloud. In a /x"*" cloud of infinite size 
equating these two indeed would seem to be the only thing that could give a 
finite result. However, one should realize that for any cloud of particles of finite 
size the standard definition of the cross section involves and integral over time 
and then the quotient of the two just mentioned times will appear in the result. 
This will depend on beam shapes. 

Secondly (2), also our result (i.e., equation 26), to be derived in the next 
section, contains a factor d^r. One could also pull a factor dt out of this and 
obtain exactly the same result as by Ginzburgs method. So, whether or not one 
likes the equating of dW/dVs mentioned in the foregoing point, one does not 
need complex momentum conservation to obtain the result that belongs to it. 

Thirdly (3), the matrix element, as seen in position space, used in Ginzburgs 
method, is not as advertised in [6]. To see this, we should realize that the method 
of Ginzburg only modifies the propagator of the muon neutrino. In particular, 
nothing is altered in the prescription of removing external propagators. The 
consequence of this is that the number of muons does not decrease. After all, 
this prescription was designed to describe a stable particle that comes in from 
infinity. This means that the production rate of neutrino's inside the ix^ cloud 
is given (in the rest frame of the ii~ ) by 

Thus, every muon produces a large number of neutrino's. Because of the con- 
servation of complex momenta, the muon neutrino gets a decay time equal to 
the decay time one would expect for the muon. The consequence is that the 
number of muon neutrino's at a particular time is given by 

N.^{t) = J dt'6{t-t')e-^^^'-''A=N,. (20) 

So in the method of Ginzburg the number of neutrino's is equal to the number of 
muons not because every muon decays into a single muon neutrino but because 
the decay constant of the muon neutrino is artificially made equal to its produc- 
tion constant. That this gives the same result as equating the two dW/dt's as 
discussed under point (1) can be easily understood because the muon neutrino's 
now only exist near the path of the fjT , so the two dfs now both refer to the 
time the jjT spends in the jjL^ cloud. 
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3.2 The Right Way 

Here we take the production process of the unstable particle into account. This 
is reasonable because unstable particles cannot really be asymptotic states. In 
the context of scalar fields, Veltman (cf., [7]) has shown that if one takes only 
stable particles for asymptotic states, one gcits a unitary S-matrix. Thus we 
must describe the unstable particle as an internal line of a larger Feynman 
graph that includes the production process. Consequently its momentum can 
have four independent components. Of course, we really do not want to include 
much information about the production process of the unstable particle in our 
calculations. For this reason we define the wave function of the unstable particle 
to be given by 



V'(P) = Y J dPadpb(l)a{Pa)4'b{pb){'2'!r)^S^{pa+Pb-'p)M{pa,Pb,p), (21) 

where A4 is the matrix element (or perhaps the sum of matrix elements) that 
describe the production process. This definition assumes that the unstable par- 
ticle in produced in a two-to-one collision, however we can safely include terms 
where the production process has as many in- or outgoing particles as we would 
like. The outgoing particles are represented by complex conjugates of the wave 
functions that they are measured to be in. Such a measurement should be car- 
ried out in order to make sure that the momentum of the unstable particle is 
fixed very accurately. When doing a calculation with this we still have to include 
the propagator of the unstable particle. The just given wave function only re- 
places the part of the matrix element that describes the production process, not 
the propagation of the unstable particle. For this propagator, the replacement 

i . J_ (22) 

— -|- imT mV 

should be used, since the momentum of the unstable particle is assumed to 
be fixed — by the just-described measurement procedure — on its mass shell so 
accurately that a function of which the width is of the order T (as the propagator 
is) cannot notice the difference. 

The normalization of this wave function, as occurs in its definition above, was 
obtained by considering the production process followed by the decay process. 
The number of events in this is demanded to be unity because we choose our 
wave function to describe one unstable particle. 

In this case the definition of the quantum distribution function should be 
modified a bit. It becomes 

n{p, r) = j ^e-^-V(p + l^^P - l^T- (23) 

Probability densities in momentum and position space can be obtained from 
this. They are given by 

/f d^P 
d^rn{p,r); p{r) = J j^--^n{p,r). (24) 

The luminosity is now dt^fined by 

dL{pi,p2,p)= ^ TO^^^r^^^'' ^l^'^^^ / d'^"rnx{p\,r)n-2.{p2,r -\- p). (25) 
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This is for the case of a stable particle colliding with an unstable one. The reader 
will presumably have no difSculty figuring out what to use for two unstable 
particles if he compares this to equation 5. The cross section is defined in 
exactly the same way as in section 2. Furthermore, the just given definition for 
the luminosity was chosen in such a way that equation 7 is kept the same too. 

After having introduced this, the calculation of the number of events pro- 
ceeds along exactly the same lines as in [2] and we will not present it here. The 
result is 

^= — / TTm'^P^ ^'^'^ / dani{pi,r)n2{p2,r + ak)fs{k\pi) 

m2 J (27r)* Jo (26) 

^J{k-P2)^ -k^pl(7{k,p2). 

The interpretation of the above formula is that the unstable particle decays at 
position r, the decay product has momentum k and then collides with the other 
particle at position r + ak, where a is a positive number, as one would expect. 
If we specialize to the case where the longitudinal beam size is much larger 
than the transversal one, we find back formula 8. The densities that occur here 
should be considered to be densities of decay events. It is perfectly natural that 
decay events are characterized by four co-ordinates. 



4 Pancake- shaped Beams 

In this section we will show that it is possible to transfer some of the imT 
to the muon neutrino propagator while keeping a clear conscience. This is in 
the context op pancake-shaped beams. In this case, we assume the transversal 
beam size to be much larger than the decay length l/F, while the longitudinal 
beam size is much smaller than the decay length. We do this by using momenta 
that have a large width in the longitudinal direction and a small width in the 
transversal direction. We assume that the unstable particle is produced at a 
position away from the position where it collides. We do this by translating the 
particle having momentum p2 in equation 1 an invariant distance s (as seen by 
the y»i-particle) away from the origin. This is done by a adding factor e^^^i'^s/™! 
in the matrix element. The complex conjugated matrix element gets a factor 
g-ispi pa/™! The variables P2 and ^i'g integration variables in momentum 
space. They integrate over the momentum distribution. 

Furthermore, we take into account (as we have done for the case of infinitely 
wide beams) the fact that unstable particles cannot really be in/out states. The 
real matrix element describing the collision process includes the production of 
the unstable particle. The incoming unstable particle differs from a stable one in 
that its momentum need not be on its mass shell. The momentum distribution 
is determined by its propagator ((p'l)^ — ml -l-imiF)"^, which makes sure that 
the momentum is peaked around the mass shell. 

We consider the quantity F which contains the factors of the matrix element 
that play a role in making the cross section finite. These are the factors that 
are peaked sharp enough to notice that p[ (and P2, of course) do not have a 
definite value but are peaked around the value pi . It is given by 

gisp2-Pl/™l 1 

^ " (P'l)^ -ml + imiT (ga - p'2)^ - M"^ + ie ^^"^^ 
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Wc now integrate over the value of (p'l)^. The integration path is chosen such 
that outgoing momenta are kept fixed, while the change A'^ of the incoming 
momenta is taken to be a linear combination of these momenta. Thus momenta 
are parameterized by 

{p',nt) = {p',nO)-tA'^; (28) 

where t is just a parameter and has nothing to do with the t-channel. Further- 
more we should demand that the value of (p^,)^ is kept fixed. It is not very 
difficult to think of a scatter setup where the particle with momentum p2 is 
a stable one, so we had indeed better not vary this one. These demands are 
satisfied by taking 

1 {pi ■P2)P2 -pIPi 

2 ipi ■ P2y - pIpI 

We take a linear combination of pi and p2 because of the pancake shape. These 
momenta only have a sufficiently large spread in momentum space in the lon- 
gitudinal direction, so this direction should be chosen for the integration path. 
Filling this into equation 27 and using contour integration over t, we find that 
for the matrix element the integration boils down to the substitution 

J.sp2-pi/m-i-sT/2 

^--2'^^A^3M^37^^(^?-m?) ^^^^ 
where a = 2 A • q2. 

We see that this no longer causes a divergence, so we will square this and 
work out the Golden Rule to find a cross section. The square of both terms con- 
tains a Breit-Wigner of which the square can be approximated hj a, 5 function. 
However, before we do this, we should discuss the meaning of the two terms. In 
the first one we effectively first integrate over the width of the unstable particle 
and then over the width of the decay product. In the second term it is the other 
way around. The interpretation is that the width that is integrated over first 
corresponds to the smallest virtuality (this corresponds to the largest distance 
scale). Thus, the first term describes events where the unstable particle does 
not decay or decays nearby the spot where the collision happens. The second 
term describes events where the unstable particle decays long before the decay 
product collides. We interpret this as the case where the unstable particle is far 
off shell and is never really produced but already decays during the production 
process. We decide to drop the second term and keep the first one. If we square 
the term that we decided to keep, we get the familiar decay law (i.e., the e~^''). 

The consequence of the just described procedure is that our squared matrix 
element contains a factor 
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For the cross section, wc find 



(5(fc' - M2). (31) 



<y{pi,P2) = f dk^^^^^^'-^a{p2,k)Uk\Pi)- (32) 



^ {pi ■ P2Y - pIpI 
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The quotient of flux factors is the gcnerahzation of the 1 — cos 9 in equation 14. 
The only factor that we might not have expected is the 1/a. 

We can shed light on this factor 1/a in the rest frame of the stable particle 
(the one with momentum p2). In that frame we have, writing out all inner 
products, 



We see that the flux factor together with the 1/a turn into a factor 1/ cos9. 
This is because we are describing the collision of two pancake-shaped beams. 
The p2 particle is taken to be in its rest frame here, so it should be imagined 
to be a pancake that docs not move. The decay product that collides with a 
particle in the pancake emerges under an angle 0. Because of this, it sees the 
thickness of the pancake enlarged by a factor 1/ cos^. 

5 Conclusions 

The problem that the cross section diverges for colliding unstable particles can 
be considered for the case of beams miich smaller than the decay length or for 
beams much larger than the decay length. We also had a look at what happens 
if two pancake-shaped beams collide. 

In the first case, the divergent part of the cross section is proportional to the 
transverse beam size. We can confirm the result of [2] that for ^/s ~ 100 GeV the 
part of the cross section independent of the beam size and the part proportional 
to it are of the same order of magnitude. However, for higher energies {^/s > 
150 GeV already) the part proportional to the beam size is at least three orders 
of magnitude smaller. At still higher energies this becomes even more, so for 
high energy colliders the linear beam size effect can safely be ignored. In practice 
this can be done by imposing cuts as one would have in a collider. 

In the second case, it should not be "solved" by introducing complex momen- 
tum components. This artifically introduces a decay time for the decay product, 
which does not seem to be a real physical effect. The authors of this paper are 
not aware of the existence of a law of "conservation of decay width" . The result 
that is obtained by the introduction of complex momentum co-ordinates can (if 
we allow for a not unreasonable modification in the definition of the cross sec- 
tion) also be obtained by using our methods. Our methods involve considering 
the wave function of the unstable particle to be a function of four momentum 
components. This is reasonable because the unstable particle cannot be an 
in/out state and really is an internal line of a bigger Feynman diagram. After 
this, the reasoning proceeds along the same lines as for the case of realistic 
beams. We find (if we do not allow for the just mentioned modification and 
assume cylindrically shaped beams) the same linear beam size effect. 

Another method, namely integrating over the width of of the unstable 
particle, is the right thing to do for the case of pancake-shaped wave packets. 

One might ask to what extent delicate gauge cancellations are destroyed 
by our approach. We do not expect these problems to be qualitatively worse 
than those encountered in, say, calculating loop corrections to LEP-2 processes 
(cf. [8]). This point will be addressed in a forthcoming publication. 




(33) 
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